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ABSTRACT: 

1. A method of measuring a flow comprising particles, 
in particular for 

blood flow measurement, in accordance with the ultrasonic 
pulse Doppler method, 

wherein ultrasonic transmitting pulses are irradiated into 
a signal entry 
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region of an object under examination in the direction of 
the particles of the 

flowing medium, and reflected echo signals received from 
the particles are 

converted into electrical Doppler signals, wherein a 
measuring signal for the 

flow of the flowing medium is determined from the frequency 
difference between 

the transmitted and received ultrasonic pulses after a 
spectral analysis, and 

wherein the measuring signal is displayed optically and/or 
acoustically, 

characterised in that it involves the following steps : a) 
a simultaneous 

spectral analysis of the Doppler signals is carried out in 
at least two 

different frequency position (alpha, beta) , where the 
Doppler spectra present 

in these frequency positions (alpha, beta) are subjected to 
laws of similarity, 

and where the spectral sector powers (MO alpha , MO beta 
and MObl , M0b2 , MOa 

) which emanate from the particles in the signal entry 
region, which move in 

the acoustic direction by means of a specific velocity 
component, are supposed 

to be identical, b) superimposed components of the Doppler 
spectra (MOa ) which 

are not to be displayed, are eliminated from the measurement 
by compensatory 

superimposition, and c) the measuring signal obtained after 
the compensatory 

superimposition is optically and/or acoustically displayed 
or recorded. 
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@ Verfahren und Oerit zur Strdmungsmessimg risen der UltraschalkPuis-Doppler-Methode. 
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@ Die Erfindung bezieht sfch auf ein Verfahren und eln 
Gerftt dteser Art, bet dem Uttraschafl-Sendepulse in ein 
Unterauchungsobjekt In Rlchtung auf eln strdmendes Me- 
dium abgestrahrt und die am strdmenden Medium reflek- 
tlerten Pulse als Echopulse empfangen werden. Aus der 
Frequenzverwerfung zwfschen ausgesendeten und empfan- 
genen Uftraschallpulsen wind eln Dopplersignal als ein Mafl 
fur die Strdmung des Mediums ermftteit Zlel der Erfindung 
1st es, ein Verfahren und ein Gerat dieser Art anzugeben, 
mrt dem slch ohne Abgang von der Puls-Doppler-Technik 
dl Berelchsgrenze der Strdmungsmessung erhdhen IdBt 
Dieses Ziel wird erflndungsgemaB wie folgt erreicht Es 
wird elne gleichzeltlge SpektraJanaryse der anfatlenden 
Dopplersignale In wenlgstens zwei unterschiedllchen, aber 
Shnlichen Frequenztagen (a, p) vorgenommen an solchen 
Stellen, an denen die Ausschnrttsleistungen (Moo, bzw. 
Mom, Moot, MoJ glelch sind; durch kompensaiorische Ober- 
lagerung werden solche Oberlagerten Arttelle der Doppler- 
spektren (M«J t die nlcht angezelgt Oder al Oberlagerter 
Bestandtell separat angezelgt werden sollen, aus der Mes- 
sung ellmlniert; das MeBergebnls wird einem Artzeigegerat 
(32, 33; 34) zur optischen und/oder akustischen Anzeige 
bzw. Aufeeichnung zugel itet 
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Verfahren und Gerat zur Stromungsmessung nach der Ultra* 
schall-Puls-Doppler-Methode 



Die Erfindung bezieht sich auf ein Verfahren zur Stro- 
mungsmessung nach der Ultraschall-Puls-Doppler-Methode 
10 gemaB Oberbegriff des Patentanspruchs 1 sowie auf ein 
Gerat zur Durchfiihrung eines sol chen Verfahrens. 

Verfahren und Gerate dieser Art gewinnen insbesondere 
fiir die Kardiologie immer mehr Bedeutung. Hier soil aus 

15 einem kleinen Signal einzugsgebiet in wahlbarer Korper- 

tiefe die Blutstromung detektiert werden. Die Brauchbar- 
keit der Ultras chall-Doppler-Technik fur die Kardiolo- 
gie ist allerdings immer noch recht eingeschrankt; aus 
energetischen Griinden konnen namlich die Parameter (z.B. 

20 Sendefrequenz u.a.) nicht so gewahlt werden, daB im vor- 
gegebenen Tief enbereich auch noch hochste Geschwindig- 
keiten der Stromung noch eindeutig erfaBt werden konnen. 
ErfahrungsgemSB f ehlt etwa der Faktor zwei fiir die Be- 
reichsgrenze der Stromungsgeschwindigkeit. Es sind zv;ar 

25 bereits theoretische Losungen des Problems bekannt; die- 
se Losungen arbeiten mit kodierten Dauerschall signal en 
in Verbindung mit einer Korrelationstechnik. Dies er- 
fordert aber das Abgehen von der bekannt en Pulstechnik, 
was nachteilig ist, weil sich der technische Aufwand er- 

30 hoht und weil eine Verbindung der Doppler-Technik mit 
der B-Scan-Technik erschwert wird. 

Aufgabe vorliegender Erfindung ist es, ein Verfahren und 
ein Gerat der eingangs genannten Art anzugeben, mit dem 
35 sich ohne Ab gang von der Puls-Doppler-Technik die Be- 
reichsgrenze der Stromungsmessung erhohen laBt. 

Kue 5 Kof / 12.08.1981 
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Die Aufgabe wird mit einem Verfahren gel 6s t, das die 
im kennzeichnenden Teil des Anspruchs 1 auf gezahlten 
Verfahrensschritte umfaBt. 

Ein Gerat zur Durchfiihrung dieses Verfahrens ist durch 
die Merkmale des Anspruchs 13 gekennzeichnet. 

Die Erfindung geht von der folgenden Erkenntnis aus: 
Dopplerspektren ein und desselben Grundgeschehens 
(Stromung im Signaleinzugsbereich) finden sich in ei- 
nem Pulssystem, das ein kammartiges Spektrum aufweist, 
vielfach in ahnlicher Form nebeneinander . Bei Bereichs- 
iiberschrei~tung des konventionellen Systems uberlagern 
sich die Dopplerspektren von Vor- und Riickstromung oder, 
falls keine Gegenstr5mung vorliegt, die Teilspektren 
der verbleibenden Stromungskomponen"ten , so daB ein kaxam 
auswertbares Gemisch entsteht. Die Ub er la gerung oder 
Uberlappung der einzelnen Spektren fallt aber in unter- 
schiedlichen Frequenzlagen des Kammspelctruins unter- 
schiedlich aus, Diese Tatsache nutz~t die Erfindung aus. 
Wird also, wie die Erfindung vorschlagt, eine gleich- 
zeitige Spektralanalyse der uberlagerten Dopplersignale 
in zwei Frequenzlagen vorgenommen, so ergibt sich ein 
MeBergebnis, das mittels Konrpensationsprinzip zur Tren- 
nung ausgenutzt werden kann, sofern die Ahnlichkeits- 
gesetze beachtet werden. Die Erfindung fuhrt also zu 
einer Trennung von Hin- und RuckfluB bei noch vertret- 
baren Fehlern durch endliche MeBzeiten. Die Bereichs- 
grenze der Stromungsmessung wird dabei betrachtlich er- 
hoht (etwa um den Faktor 5/3 = 1,67), wobei von der 
Pulstechnik nicht abgegangen zu werden braucht. 

Veitere Vorteile und Einzelherten der Erfindurg c-rgebcn 
sich aus der nachf olgenden Beschreibung eines A-?sfUh- 
rungsbeispiels anhand der Zeichnung und in VerV uiung 
mit den Unteranspruchen . ^ ORIGINAL 
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Es zeigen: 

Figur 1 ein Diagramm, das fur zwei unterschied- 
liche Frequenzlagen die Xhnlichkeitsge- 
setze zwischen den Leistungsspektren 
einer StrSmungsrichtung veranschaulicht , 

Figur 2 das Prinzip einer kompensatorischen Mehr- 
bereichs-Frequenzanalyse fiir iiberlagerten 
Hin- und RuckfluB , 

Figur 3 das Prinzip der Figur 2, angewendet auf 
die Randbereiche der Dopplerspektren, 

Figur 4 das Prinzip der Aufaddition des isolier- 
ten Spektrums, 

Figur 5 ein Ausfuhrungsbeispiel fiir eine kompen- 
satorische Mebrbereichs-Frequenzanalyse 
im Frinzipschaltbild, 

Figur 6 einen Signal-Synthetisierer zur Erzeu- 
gung eines kiinstlichen Doppler signals 
mit der Leistungsdichtefunktion eritspre- 
chend dem Frequenzanalyseergebnis bei 
iiberlagerten Spektren. 

In der Figur 1 sind zwei Dopplerspektren in zwei unter- 
schiedlichen Frequenzlagen und (i ausschnittsvergros- 
sert dargestellt. Die unterschiedliche Lage im Spektrum 
entspricht der gleichzeitigen Untersuchung des Slgnal- 
einzugsgebiets mit verschiedenen Send ekreisfr ecu unzen 
SL^ und il y Die gepunkteten Rechteckflachen rc-pra- 
sentieren gleiche Leistungen M Qc< und M Q ~ . Sie sii:d 
in den Mittenfrequenzen Ct ^ bzw. K • M ^ und .in ihren 
Bandbreiten j^Cj bzw. K • A CO entsprechend dem ^oh'-vo^s- 
faktor K = verandert. Dasselbe gilt fiir die ) t Oi -v.er- 

'BW) ORIGINS 
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ten Mittenfrequenzen — — — - bzw. — = — — und fiir die 

.Si -j jZ. 5 

normierten Bandbreiten bzw. K - A O ^ jgine umge- 

kehrt proportional e Beziehung ergibt sich zwischen den 

2 

Leistungsdichten A bzw. g- der Leistungsspektren, so 
daB die Leistungs- bzw. Flachengleichheit M Q ^ = M Q fi 
eingehalten wird. Damit sind durch Figur 1 die fiir bei- 
de Leistungsspektren giiltigen Ahnlichkeitsgesetze de- 
finiert. Ahnliche Leistungsspektren haben jedoch nicht 
gleicbzeitig auch identische Zeitfunktionen f-j(t) bzw/ 
f^Ct) des Dopplersignals zur Folge. Diese sind im all- 
gemeinen unterschiedlich. 

Die Figur 2 zeigt in Anwendung der in der Figur 1 ver- 

anschaulichten Ahnlichkeitsgesetze das Prinzip einer 

kompensatorischen Mehrbereichs-Frequenzanalyse fiir 

uberlagerten Hin- und RiickfluB. In der Figur 2 sind 

im Diagramm oben jeweils wieder in den beiden Frequenz- 

lagen oc und |S je zwei sich uberlagernde Spektren dar- 

gestellt. Das Spektrum | F^ ^ } aus der Lage c* und das 

Spektrum I ^ I ^ aus der Lage £ stammen dabei vom Hin- 

fluB. Die Spektren |F atf | 2 der Lage c* und lF a p,) 2 der 

Lage £ result! eren entsprechend aus dem RiickfluB. Mit- 

tels zwei FFI-Spektralanalysatoren FFT1 und FFT2 mit 

den Filterfunktionen und Gp> der variablen Mitten- 

frequenz Cc^ 9 K> und DurchlaBbandbreiten A CJ und 

K • A CD werden in den beiden Frequenzlagen und ^ 

aus den Spektren des Hin- und Ruckflusses unter Beriick- 

sichtigung der Ahnlichkeitsgesetze die Summenleistungen 

bzw. Ausschnittsmomente 0. Ordnung herausgefiltert. Im 

Falle d es Spektralanalysators FFT1 erhalt man die Aus— 

schnittsleistung £ Ifl = M , 0 + M . Der FFT2 llefert 

ooc ooZ oa 

die Ausschnittsleistung IM = M , „ + M . Nach don 

° o odi oa 

vorstehend abgehandelten Ahnlichkeitsgesetzen sl^d die 
Leistungsausschnitte M Qb1 , M Qb2 und M Qa fiir beidc Fre- 
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quenzlagen oC und £ gleich, was im Diagramm in der Fi- 
gur 2 oben dargestellt ist. Werden also die Ausgangs- 
signale der beiden Spektralanalysatoren FFT1 und FFT2 
voneinander subtrahiert (in einem Subtrahierglied D) , 
so erhSlt man als Ergebnis die Leistungsdiff erenz 
M ol;)1 - M Ql3 2 aus dem Spektrum des Hinflusses. Der Lei- 
stungsanteil aus dem RuckfluB ist hingegen eliminiert. 
Die Anteile der erhaltenen Leistungsdiff erenz liegen, 

bezogen auf (& 09 im Frequenzabstand = 2 ^T • 

K-1 * 
-g- auseinander, wie in der Figur 2 unten angedeutet 

ist. 



Die kompensatorische Mehrbereichs-Frequenzanalyse er- 
laubt also die Unterdruckung einer FluBrichtung. Das 
1 5 Spektrum der nicht unter driickt en FluBrichtung, d.h. 
das isolierte Spektrum, erhalt man in einfacher Weise 
unter Anwendung des Prinzips der Figur 4 durch Auf- 
addition nacheinander in Schritten I, II etc. erhal- 
tener Werte, beginnend mit dem Schritt I vom Rand der 
Spektren, bei sich schrittweise anderndem Co ^ bzw. 
K*0j 1 der Spektralanalysatoren FFT1 und FFT2. 



Die Randwertanalyse ist in der Figur 3 dargestellt. In 
diesem Falle lief ert der Spektralanalysator FFT1 im 

25 Schritt I an seinem Ausgang lediglich einen Leistungs- 
anteil M Qa , da Ja das Spektrum If^^ 2 in diesem Fre- 
quenzrandbereich bei der fiir den Hinflufl maximal er- 
laubten Bandbreite Null ist. Der Spektralanalysator 
FFT2 lief ert hingegen im Randbereich in der Frequenz- 

30 lage K co-j unter Beriicksichtigung der ftinlichkeitsge- 
setze das Ergebnis M Qb1 + M Qa . Nach Diff erenzbildung 
im Diff erenzbildner D ergibt sich die Leistung M Qt)1 . 
Beginnend von dieser Randleistung des Randanalysier- 
schrittes I konnen nun entsprechend dem Diagramm der 

35 Figur 4 die sich ergebende Diff erenzleistung M ob1 - 
M ob2 der Fi S ur 2 bzw - bei Durchfiihrung entsprechend 

BAD ORIGINAL 
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weiterer Analysenschritte- II, III, XV etc. weitere Dif- 
ferenzleistungen aufaddiert werden. Es ergibt sich dann 
das in der Figur 4 dargestellte isolierte Leistungs- 
spektrum |F b | 2 des Hinflusses. Das Spektrum des Ruck- 
5 flusses|F a J 2 - kann in ahnlicher Weise, beginnend vom 
linksseitigen Randwert, und durch entsprechende Diff e- 
renzbildung mit den Leistungswerten des HinfluBspektrums 
gemessen werden. 

10 Der vorstehend beschriebene Sachverhalt wird im nach- 
folgenden durch eine mathematische Leistungsbetrachtung 
noch mehr vertieft: 

Entsprechend den Figur en 2 und 3 konnen in den unter- 
15 schiedlichen Frequenzlagen c* und (S also durch einen 
Einseitenbandempfanger die folgenden Spektren empfan- 
gen werden: 

niedrige Frequenzlage (F &oC ( c^) + F b ( 0) 2 ) )(1) 

20 - - • 

hohere Frequenzlage < F a F> ( ^ 3> + F bf» ( W 4> 

Diese uberlagerten Spektren werden nun einer Filterung 

in Kombination mit einer Lei stungsmes sung unterworf en. 

25 Die Filterfunktionen G ^ und G ^ sind so gewanlt, daB 

gemaB den vorausgesetzten Shnlichkeitsbedingungen sich 

entsprechende Spektralanteile in F und F „ elei- 

a oc a 6 6 

chennaBen erfaBt werden. Man stelle sich also vor, durch 
Einseitenbandtechnik liegen die beiden Frequenzlagen or 

30 und £ z.B. mit Bezug auf Sl^ bzw. SL^ demoduliert 

gleichzeitig vor. Zwei Spektralanalysatoren bilden die 
Leistungen aus der Bandbreite ^ « bei der Frequenz C?^ 
hzw. K • A co bei der Frequenz K • Ci v Man erhalt dann 
die Leistungen bzw. Ausschnittsmomente 0. Ordnung mit 

35 den Empfangsempfindlichkeiten « und p fur diepe Be- 
reiche: 
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£ M 



oet 

o 



"J J F aw < «« 1> # G «*I + 



2JL 2 



(3) 



10 

2^ 



15 



o 

2 W 



(4) 



"A 



20 Mi-t lG ec I s 1 bzw. lG p \ = 1 im DurchlaBbereich ergibt 



sich dann l>ei kleinem A OJ und |F av (CD-) " G of | - 



2 5 G ( i l 2 = A f> 2 » K p (C i 4 , - G p| 2 - 2 



; (o/ a+ * B «; (cc 2 ) Ato (5) 

30 bzw. mit 

2 

A p;(KOi 1 ) ^ (6) 
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Aus den Gleichungen 5 und 7 ergibt sich dann mit folgen- 
der Zwischenrechnung 



(8) 



wegen 



B % (K O 2 ) - £ b2 « ( 2 ) (9) 



i M o« -f^op =^AOK (B^ ;(K052) 



b 2 ^ . _ K cj .,) ) - Aca K oc A b| (10) 



oder 



AUc<Ab5 

Diese Gleichungen sagen aus, daB die Differ enz der Lei- 
s*ttmgen gemessen werden kann, welch e in den Eir-zel spek- 
tren F^ ^ bzw. ^ an den Stellen O 2 hzw. K Oi 2 

(^=ST" - K 0i 1 ) bzw. (f^r- - 63 1 ) mi-t der Bandhreite Aca 
hzw. K AO ausgefiltert werden. Die Differenz der Fil- 
"termittenfrequenzen ist in der Frequenzlage : 



(12) 
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Bei einer zentralen Sendefrequenz von 3,7 MHz, einer 
Pulswiederholfrequenz von ^ » 15 kHz und einem Abstand 
zwischen den Frequenzlagen e< und £ von 660 kHz er- 
gibt sich dann z.B. 

5 

Af K£ = 2 ' 97 * 3 K ^ * (14) 

Bezogen auf das Frequenzband ex. gilt 

10 Af K<* - A K K * - 2 » 4 79 « 2,5 KHt (15) 

Man kann dies als eine Frequenzanalyse von F b ^ mit der 
Bandbreite Ao auffassen, bei der die Amplituden des 
Leistungsspektrums an zwei Stellen im Spektrum heraus- 
15 gegriffen und die Differenz gebildet wurde. 

■ 

Die Filtermittenfrequenzen beider Stellen haben konstan- 
ten Abstand, der von den Systemparametern *t? und K ab- 
nMngt (Gleichungen 12 und 13). Die Filtermittenf requen- 
20 zen sind . 

furF bp : oben: K Ci 2 = K C 2 -^— C0 1 ) (16) 



25 



30 



unten: 



2TT v 2 
X 



"KCi 1 = - ■=-£ — (K -1) + 

1 C 



KU 2 = KG> 2 - Ato Kp> 



(17) 



fur F b ^ : oben: c* z = IX- - Co 1 ( 18) 
unten:. - « ^ _2X.Ll! + C02 = 



M . (19) 
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Die Differenzfrequenz Af Ko ^ ist genau so groB wie di 
Frequenzliicke zwischen der maximalen Frequenz von F bcC 
und 1/-£ f die bestehen bleiben muB, damit ^ nicht 
iiber die Bandbreite wachst (siehe Figuren 2 und 3): 

5 . 

Wenn in der Frequenzlage [£ gilt 

K °2 - ^ (20) 
10 dann ist die Frequenzliicke in der Lage c< 

(21) 

15 Diese Tatsache kann man nut z en, urn aus den Dif f erenzwer- 
ten der Spektralanalyse auf das Spektrum selbst zu 
schlieBen: 

Legt maxt die obere Filterstelle, bezogen auf F b( ^ , in ' 

20 den Bereich G) * = 0 bis Af 7rf bzw. C3 Q = 4=2 A f-^ 

bis dann *ird |F b<x - G^ 2 = 0 bzw. B 2 ^ ^ ^ = 

0, wenn keine Bereichsuberschreitung vorliegt (s. Fi- 
gur 3). In diesem Bereich kann also B^. ^T ^ ^ di- 
rekt gemessen werden, und zwar an einer 'ffei waMl^aren 

25 Stelle mit einstellbarer Bandbreite je nach Filterwahl. 
Im Beispiel bedeutet das, bezogen auf F^ ^ , daB der 
2,5 kHz breite (hochstzulassige) Bereich von 10 kHz bis 
12,5 kHz bzvr. -bezogen auf F^ ^ der 3 kHz breite Bereich 
von 12 bis 15 kHz direkt ausgemessen werden kann. Die 

30 darunterliegenden Bereiche findet man durch Aufaddition 
der Differenzwerte (s. Figur 4). Im Beispielfall mit 
A f Kdi =2,5 kHz sind es vier Aufadditionen fiir das 
gesamte Spektrum. 
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Das Spektrum der GegenstrSmung findet man nach Glei- 
chung (5) durch eine dritte Spektralanalyse z.B. von 
M o<* 9 wobei fiir gleichartig liegende Filter G ^ gilt 

5 * A oC ; (*) 1 ) A » = * M o* (ts n ) -^(U^AO (22) 

Die Analysebreite A« ist als frequenzunabhangig kon- 
stant angenommen. 

10 Das Prinzipschaltbild zur Durchfuhrung der kompensatori- 
schen Mehrbereichs-Frequenzanalyse ist in der Figur 5 
dargestellt. In diesem Prinzipschaltbild erzeugt ein 
Oszillator 1 die zentrale Sendefrequenz SL/2rK (z.B. 
3,7 MHz) xmd Hilf sfrequenzen, wie die Pulswiederholfre- 

15 quenz 1/t (z.B. 15 kHz) und Empf angs-Oszillatorfre- 
quenzen ^/2Tf ; , sovie die Empfangstorzeit 

T £ (z.B. 2,7/us) und Zeitverzogerung ■£ 2 fiir die Tiefen- 
lage des Signaleinzugsbereichs. In einem Sendetor 2 
wird eine Zahl von Perioden der Sendefrequenz (z.B. 

20 1 Periode) im Rhythmus 1/*T ausgeblendet uiid iiber den 
Sendeverstarker 3 zum elektroakustischen Wandler W ge- 
schickt. 

Das Empfangssignal gelangt iiber einen Empfangsverstar- 
25 ker k zu zwei Einseitenbandempfangern 5 und 6. Diese 
empfangen ein Frequenzband der Breite 1/t in einer 
tiefen Frequenzlage ©t (z.B. unterhalb -^-,) und einer 
hohen Frequenzlage p» (unterhalb Die Einseiten- 

bandempf anger konnen z.B. nach der Quadra turmethode ar- 
30 beiten. Die Ausgangssignale der Empf anger sind die de- 
modulierten Dopplersignale, welche von Vor- und Ruck- 
fluB herruhrende uberlagerte Spektren beinhalten konren. 

Von beiden Aus gangs signal en wird in den Integra tor en 7, 
35 8, 9, 10 (Serienschaltung aus Quadriergliedern und Tief- 
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passen) die jeweilige Gesamtleistung (Zeitkonstante z.B. 
1/5 Hz = 0,2 s) und im Quotient enbildner 11 der Quo- 
tient */ f beider Werte gebildet. 

r 

5 Beide Ausgangssignale werden auBerdem fiinf FFT-Spektral- 
analysatoren FFT1 bis FFT5 (BIScke 12 bis 16) angeboten. 
Die Analysebandbreite fiir das Signal der Frequenzlageoc 
ist z.B. A.f = 125 Hz, FUr die Frequenzlage g> betragt 
sieKAf = 150 Hz, wobei = ^3/ SL -, = 1,2. 

10 

Zuerst soil das Zusammenwirken von FFT2 und FFT4 be- 
schrieben werden: Wahrend FFT2 die spektrale Leistung 
bei der Kreisfrequenz CJ 1 bildet, berechnet FFT4 die- 
jenige bei Kto^ Das Ergebnis von FFT2 wird in einem 

15 Subtrah i er glied 19 von dem im Multiplizierglied 17 mit 
**/ (S multiplizierten Ergebnis von FFT4 abgezogen. Auf 
diese ¥eise verschwindet der von | F & | 2 herruhrende Teil 
des Spektrums und es bleibt die Differenz zweier Spek- 
tralteile von j F^J , die im Abstand A ± K liegen. Solange 

20 0< CJ^ «<f A ^ Koc * wird das Leislningsspektrum "von F^ 
unmittelbar im Bereich I (s. Figur 4) gebildet, wenn 
keine unzulassigen Bereichsttberschreitungen der Band- 
breiten bzw. Ge s chwindi gkei t en vorliegen. Venn 
A.CJ K ^<" 1 <2AO Kc< , wird der Bereich II 

25 von F^ erfaBt und man kommt zum Analyse ergebnis durch 
Aufaddition (Addierglied 24 mit vorgeschaltetem Spei- 
cher 22) der Werte mit dem vorher fiir Bereich I erhal- 
tenen Wert. Gabe es keine (meBzeitbedingten) MeBfehler, 
konnte in dieser Weise fortgefahren werden, bis das 

30 ganze Spektrum von F b ^ analysiert ist. Kleinere (noch 
tolerable) Fehler erhalt man, wenn die (noch f ehlenden) 
Bereiche III, IV, V von F^ auf andere ¥eise berechnet 
werden. Die Analysatoren FFT2 und FFT4 werden also nur 
uber die Bandbreiten 2 A * K(7 < (2 • 2,5 kHz) bzw.2Af R ft 

35 (2-3 kHz) benotigt. 
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Entsprechend arbeiten die Analysatoren FFT3 und FFT5 
zusammen. Sie arbeiten bei den Kreisfrequenzen u) 2 bzw. 
k 2 auch nur Uber die Bandbreiten 2 A f R bzw. 
2Af K(i < von %f - 2 A f Kef . ^ bis 2|L) ^ bilden 
5 nach Multiplikation im Multiplizierglied 18 und Subtrak- 
tion im Subtrahierglied 20 sowie Auf addition der Diffe- 
renzwerte im Addierglied 25 das Spektrum von F & ^ in den 
Bereichen IV, V. Das Spektrum. von F b<< in diesen Berei- 
chen IV, V findet man durch Diff erenzbildung im Diffe- 

10 renzbildner 26 mit dem Analyseergebnis -von FFT1 (Kreis- 
frequenz Bandbreite 0 bis 2 ^-), welches die tiber- 

lagerten Signale betrifft. Der unterste Bereich V von 
F b < ist Sleich dem Ergebnis von FFT1 , weil F & ^ nicht 
in diesen Bereich vordringt, wenn keine unzulassige -Be- 

15 reichstiberschreitung vorliegt. Dasselbe passiert im 

Differenzbildner 27, um die noch fehlenden Bereiche von 
F a oc zu f inden. Vom Auf wand her entsprechen FFT4 und 
FFT5 zusammengenommen einem Analysator der Art von FFT1 . 
Die Ergebnisse von FFT2 und FFT3*k6nnen auch direkt von 

20 FFT1 erhalten werden. . Somit sind also vom Gesamtauf wand 
her gesehen zwei FFT -Analysatoren nach Art des FFT1 
hardwaremaflig notwendig. Die Aufteilung auf fiinf FFT- 
Prozessoren mit dem Gesamtaufwand von eigentlich nur 
zwei FFT-Prozessoren wurde aus GrUnden der besseren 

25 Obersichtlichkeit gewahlt. Auch die BlScke 7, 8, 9, 
10 miissen nicht hardwaremSBig vorhanden sein; sie 
kSnnen durch Summenbildung der Analyseergebnisse bei- 
der Spektrographen ersetzt werden. 

30 Im Block 28 wird ein Naherungswert fur die Integral- 
funktionen 



J, 



bor /l F b*l 2dw 2 = i U^lMo, (23) 



35 co 2 
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(24) 



durch Bildung der Summen der Analyseergebnisse gebildet 
und anhand der Funktionen entschieden, ob spektrale 
Uberlappung beziiglicli der Frequenzlage «r vorliegt. Wenn 
dies nicht der Fall ist, wird das Analyse ergebnis von 
FFT1 direkt verwendet, urn die geringstmoglichen Fehler 
zu erhalten. 

Anhand des Endwertes der Integrale ( Q = O 2 = /t> ) 
kann auch entschieden werden, ob ein Spektrum gegeniiber 
den anderen so schwa ch ist, daB mit unzulassig groBen 
Fehlern gerechnet werden muB. Man kormte dies im Anzei- 
gegerat 34 anzeigen, indem zJB. die eingeblendete Kur- 
ve f bei zu schwa ch em Spektrum weggelassen wird. 



In Mittelwertbildnern 29, 30 werden mehrere benachbarte 
Analysewerte im Sinne einer Mittelung zusammengefaflt, um 
zu weniger stark streuenden ¥erten zu kommen. Die Zahl » 
der zu mittelnden Werte richtet sich nach der Frequenz- 
auflosung, die. z.B. 10 % der Frequenz betragt. j 

Damit stehen die Einz el sp ektr en trotz gegensei tiger tlber- 
lagerung getrennt zur Weiterverarbeitung zur Verfugung. 
Es bietet sich eine optische und akustische Anzeige an. 

Die akustische Anzeige ist kein Problem, wenn in Block 28 
festgestellt wurde, daB keine Uberlappung der Spektren • 
vorliegt. Dann wird in den oberen Teilen d r Blocke 32, 
33 eine Trennung der Signale uber (variable) Filter vor- 
genommen und das Band | F b< ^| durch Demodulation nsit der 
Oszillatorfrequenz invertiert; anschlieBend vc-rden 

beide Richtungs signal e getrennt wiedergegeben. 

- BAD ORIGINAL 
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Im Falle uberlagerter Spektren sind die Zeitfunktionen 
zu den Einzelspektren F &c< und ^ aus | F & ^ | 2 bzw. 
l F b*.l nic bt rekonstruierbar. BekanntermaBen reprasen- 
tieren aber nur "grobe Ausschnitte" des Dopplerspektrums 
gemaB der Frequenzauflosung des Puls-Doppler-Systems ei- 
ne Aussage ttber die Stromung. Die Doppler-Zeitfunktionen 
ein und desselben Stromungsfadens, aus verschiedenen 
Frequenzlagen gewonnen, unterscheiden sich stark. 

Es muBte also erlaubt sein, ein Rauschspektrum gemaB 
I F aoi l 'bzw. lF bo( i 2 neu zu generieren und akustisch 
als Ersatzfunktion wiederzugeben, ohne dadurch einen 
wesentlichen Informationsverlust bezUglicb der Stro- 
mungsaussage hinnehmen zu muss en. Die Zeitfunktion hat 
dann mit dem ursprlinglichen Dopplersignal nur noch ihre 
"spektrale Eigenschaft" gemeinsam. 

Daher symbolisieren die unteren Teile der Blocke 32, 33 
also Signal-Synthetisierer, z.B. gemaB Figur 6, die umge- 
kehrt wie ein Analysator arbeiten. Im vorstellungsmaBig 
einfachsten Fall besteht er aus einer Filterbank 35 bis 
39 etc. (entsprechend den Analysebereichen) mit Stell- 
gliedern 40 bis 44 etc. fiir die Verstarkung jedes Fil- 
teraweiges. An alien Eingangen liegt ein breitbandiges 
weiBes Rauschen eines Rausch genera tors 45. . Die Aus gang e 
werden alle mittels Summierschaltung 46 aufsummiert. 
Die Stellglieder werden gemaB dem Analyseergebnis ein- 
gestellt. Technisch elegant ere (digitale) Auf fuhrungs- 
formen sind denkbar. 

Die optische Anzeige der Spektren kann in bekannter Wei- 
se erfolgen (mittels Anzeigegerat 34), wie insbesondere 
Kathodenstrahlrohre. Die Analysefrequenzen werden iiber 
der Zeitachse dargestellt, vobei die spektrale Lei- 
stungsdichte bei jeder Frequenz in Helligkeitsmo duration 
gebracht wird. Zusatzlich kann eine Schatzfunktion fur 



5 



10 



20 
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die mittlere Frequenz f und die Bandbreite <s* F einge- 
blendet werden. Die mittlere Frequenz erhalt man aus 
dem Einzelspektrum nach der folgenden Formel, indem man 
die Integrale durch endliche Summen ersetzt. 

oo 

* M .oo c M J ' F («)| 2 «'4« 

Co = 2 TT f = 2 2 (25) 



• lim rT ( l y («)l 2 «di 

M"*°° M o 



Die Bandbreite berechnet man ebenfalls aus dem Spek- 
trum nach der Fonnel: 



^ /Ml - 2 " 

15 F = "f fiT " f (26) 

wobei M Q - 0. Moment des Dopplerspektrums 
M2 = 2. Moment des Dopplerspektrums 



Die am Anzeigegerat 34 angedeuteten Kurven stellen 



t - 6~ F dar. 



6 Figuren 

24 Patentanspriicbe 



Pa tent an sp ruche 
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1 . Verfahren zur StrSmungsmessung nach der Ultraschall- 
Puls-Doppler-Methode, bei dem Ultraschall-Sendeimpulse 

5 in ein Untersuchungsobjekt in Richtung auf ein stromen- 
des Medium abgestrahlt und die am stromenden Medium re- 
flektierten Pulse als Echopulse empfangen werden und 
bei dem aus der Frequenzverwerfung zwischen ausgesen- 
deten und empfangenen Ultra schallpulsen ein Dopplersi- 
10 gnal als ein MaB fur die StrSmung des strSmenden Me- 
diums ermittelt wird, dadurch gekenn- 
zeichne-t, daB es die folgenden Verfahrens- 
schritte aufweist: 

15 a) es wird eine gleichzeitige Spektralanalyse der 
anfallenden Dopplersignale in wenigstens zwei 
unterschiedlichen, aber ahnlichen Prequenzlagen 
( ** t P ) vorgenommen an solchen Stellen, an 
denen die Ausschnittsleistungen (M , M bzw. 

20 M QbV , M qS2 , M oa ) gleich slnd, °* °P 

b) durch kompensatorische Uberlagerung werden sol- 
che tfberlagerten Anteile der Dopplerspektren 
( M oa )» d±e nicn "t angezeigt oder als tiberlagerter 

25 Bestandteil separat angezeigt werden sollen, aus 
der Me s sung eliminiert, 

c) das MeBergebnis wird. einem Anzeigegerat (32, 33; 
34) zur optischen und/oder akustischen Anzeige 

30 bzw. Aufzeichnung zugeleitet. 

2. Verfahren nach Anspruch 1, dadurch g e - 

k e n n z e i c h n e t , daB die Frequenzanalysen in 
den unterschiedlichen Frequenzlagen ( t* , f ) an sclchen 
35 Stellen vorgenommen werden, an denen sich Ausschniits- 
leistungen ergeben, die in den Mittenfrequenzen ( O 



fiM> ORIGINAL 



0075195 



- 18 - VPA 81 P 5076 E 

bzw. K • ta -j) und in ihren Eandbreiten (Aco bzw. K -Acj ) 
entsprechend einem Dehnungsfaktor K verandert sind und 
bei denen slch betreffend diesem Faktor K eine umge- 
kehrt proportionale Beziehung zwischen den Leistungs- 
dicht^i (A bzw* g-) der Leistungsspektren ergibt, mit 
K = — , wobei SL* undXL verschiedene Sendekreis- 
frequenien in den unterschiedlichen Frequenzlagen ( K 
und f» ) darstellen. 

3- Verfahren nach Anspruch 1 oder.2, d a d u r c h 

gekennzeich.net , daB bei gleichzeitig 

erf aBtem Hin- und Ruckf luB des stromenden Mediums ei- 
2 

nes (lF"b| ) der uberlagerten Spektren, resultierend 
aus der einen FluBrichtung, isoliert werden kann da- 
durch, daB beginnend von einem Randwert (Figur 3) , bei 
dem wenigstens der Leistungsanteil des zu isolierenden 
Spektrums in der einen Frequenzlage ( K ) Mull ist, 
nach der kompensatorischen Methode zuerst der Randwert 
(M Q ^)j des zu isolierenden .Spektrums aus der zweiten 
Frequenzlage ( £ ) gemessen wird und anschlieflend durch 
schrittweise WeiterfUhrung der Analyse (Figur 2) mit 
sictL andernder Doppler-Frequenz weitere Diff erenzwerte 
des isoliert en Spektrums an zwei unterschiedlichen 
Frequenzstellen ermittelt und sukzessive zu vorher er- 
mittelten Werten aufaddiert werden (Figur 4). 

4. Verfahren nach Anspruch 5, dadurch g e - 
ke.nnz e i chne t , daB eine Analyse, beginnend 
von einem Randwert urn eine vorgebbare Zahl von Schrit- 
ten in der einen Richtung, von einer Gegenanalyse urn 
eine entsprechend vorgebbare Zahl von Schritten, be- 
ginnend vom anderen Randwert in der Gegenrichtung r be- 
gleitet wird. 
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5. Verfahren nach Anspruch 3 oder 4, dadurch 
gekennzeichnet, daB das Spektrum | F a | 2 
aus der GegenfluBrichtung isoliert werdenkann dadurch, 
daB die Differenz zwischen dem tiberlagerten Spektrum 

5 (Ausgangssignal von FPT1 in Figur 5) und dem zuerst 
isolierten Spektrum lF b | 2 gemessen wird. 

6. Verfahren nach Anspruch 3 oder 4, dadurch 
gekennzeichnet, daB das Spektrum { F | 2 

10 aus der GegenfluBrichtung isoliert werden kann dadurch, 
daB anstelle der zuvor ausgemessenen Seitenbander die 
jeweils anderen Seitenbander in beiden Frequenzlagen 
ausgewertet werden. 

15 7. Verfahren nach einem der Ansprtiche 3 bis 6, da- 
durch gekennzeichnet, daB vor 
Durchftihrung einer Analyse von den beiden tiberlagerten 
Spektren aus Hin- und RiickfluB in den beiden unter- 
schiedlichen Frequenzlagen ( oc , f> ) die Gesamtlei stung 

20 der tiberlagerten Spektren gebildet wird (in Leistungs- ' 
bildnern 7, 9 bzw. 8, 10 gem&B Figur 5) und daB die 
Information tiber die Gesamtleistungen dazu verwendet 
wird, die Pegel der Dopplersignale, die den beiden 
Frequenzlagen entstammen, elnander anzugleichen. 

25 

8. Verfahren nach Anspruch 7, dadurch ge- 
kennzeichnet, daB der Pegelangleich in 
der Weise geschieht, daB die ermittelten Leistungen 
der tiberlagerten Gesamt spektren durcheinander dividiert 
30 werden (im Dividierglied 11 gemaB Figur 5) und an- 
schlieBend das Doppler-Signal eines Kanales mit dem 
gebildeten Quotientensignal multipliziert wird (z.B. 
in den Multipliziergliedern 17, 18 gemaB Figur 5). 

35 9. Verfahren nach einem der Anspriiche 1 bis 8, da- 
durch gekennzeichn t, dafl die 
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Analyse mit einer absoluten Bandbreite (Au, KAco) 
vorgenommen wird, die der Bandbreite entspriclrt, welche 
bei der niedrigsten interessierenden Doppler-Frequenz 
gemaB der Unscharfer elation des Doppler-Systems In- 
5 formation Uber die entsprechend interessierende 

niedrigste StrSmungsgeschwindigkeit tragt (entsprecliend 
der MeBgenauigkeit fur die Geschwindigkeit). 

10. Verfahren nach Anspruch 9, dadurch g e - 
10 kerinzeicJanet, daB filr hShere interessieren- 
de Doppler-Frequenzen solch eine Anzahl von aus ent- 
sprechenden Doppler-Frequenzbereichen stammenden 
Analy sewer ten der absoluten Bandbreite (AO , KAA>) 
aufaddiert (in den Blocken 29, 30 gemaB Pigur 5) wer- 
15 den, die in der Summe die Bandbreite ergeben, welche 
bei der h5heren interessierenden Doppler-Frequenz 
gemaB der UnschsLrf erelation des Doppler-Systems In- 
formation uber die entsprechend hohere Stromung tragt • 

20 11. Verfahren nach einem der Anspruche 1 lis 10, d a - 
durch gekennzeichnet, daB aus dem 
Analyseergebnis riickwirkend ein Doppler-Signal 
synthetisiert wird (in den Rauschsynthetisierern, 
die in den Blocken 32, 33 der Figur 5 enthalten bzw. 

25 in Figur 6 detailiert dargestellt sind), dessen Zeit- 
funktion zwar nicht der jenigen entspriclrt, die im iiber- 
lagerten Spektrum enthalten ist, die aber die.gleichen 
spektralen Eigenschaften (Leistungsspektrum) wie diese 
Zeitfunktion hat. 

30 

12. Verfahren nach einem der Ansprttche 1 bis 11, d a - 
durch gekennzeichnet, daB in eine 
helligkeitsmodulierte Spektraldarstellung des Spektrums 
(Frequenz uber der Zeit) als Kurvenzug eine mittlere 
35 Frequenz (f ) des Analyse ergebnisses gegebenenfalls 
zusammen mit Streuwerten einblendbar ist (Signalauf- 
zeichnung auf dem Auf zeichnungsgerat 34 gemaB Figur 5). 
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13. Gerat zur DurchfUhrung des Verfahrens nach einem 

d r Anspriiche 1 bis 12, dadurch gekenn- 
zeichnet , daB zur Spektralanalyse der anfal- 
lenden Doppler-Signale in wenigstens zwei unterschied- 
5 lichen, aber ahnlichen Frequenzlagen wenigstens 
zwei Spektralanalysatoren (FFT1 , FFT2 bzw. FFT1 bis 
FFT5) vorhanden sind, die Filterfunktionen (G K , Gp ) 
sowie variable Mittenfrequenzen ( CO 1 , K •ts^) und 
DurchlaBbandbreiten 9 k •a.o) derart aufweisen, 

10 daB aus den beiden unterschiedlichen Frequenzlagen % 

(<*, f> ) vor oder nach Demodulation unter Berticksichti- 
gung der JShnlichkeitsgesetze. gleichartige Summen- 
leistungen bzw. Ausschnittsmomente herausgef iltert 
werden. 

15 

14. Gerat nach Anspruch 13, dadurch g e - 
kennzeichnet, daB die Svmnnenleistxangen 
bzw. Ausschnittsmomente an den Ausgangen der Spektral- 
ana l y satoren , nach vorzugsweise vorhergehendem Pegel- 

20 ausgleich, einem Subtrahierglied (D) zur Bildung der 
Differenz zugeleitet werden. 

15. Gerat nach Anspruch 14, dadurch g e - 
kennzeichnet , daB das Diff erenzergebnis 

25 einem Anzeigegerat (32 r 33; 34) zur optischen und/oder 
akustischen Anzeige bzw. Aufzeichnung zugeleitet wird. 

16. Gerat nach einem der Anspriiche 13 bis 15, da- 
durch gekennzeichnet, daB die 

30 Spektralanalysatoren in der Weise ausgebildet sind, 
daB eine schrittweise Analyse, beginnend von einem 
Randwert, ermoglicht ist. 

17. Gerat nach Anspruch 16, dadurch g e - 

35 kennz ichnet, daB eine Aufaddierschaltung 
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(24, 25) vorhanden ist, die die Differenzwerte der 
Ausschnittsleistungen, di bei schrittweiser Weiter- 
fuhrung der Analyse erhalten werden, sukzessive zu 
zuvor ermittelten Werten aufaddiert, wodurch das 
5 Lei stung sspekt mm einer FluBrichtung erhalten wird. 

18. Gerat nach Anspruch 17, d a d u r c h ge- 
kennzeichnet, daB die Spektralanalysatoren 
auch in der Weise ausgebildet sind, daB eine schritt- 

10 weise Analyse gleichzeitig aus zwei entgegengesetzten 
Frequenzrichtungen, beginnend von jeweils gegeniiber- 
liegenden Randwerten der Spektren, mBglich ist. 

19. Gerat nach einem der Anspriiche 13 bis 18, d a - 

15 durch gekennzeich.net, daB oe einem. 
Einseitenbandempf anger (5, 6) fur zwei' Pre quenzlagen 
(-f, {!>) wenigstens je ein Spektralanalysator (FFT2, 
FFT3 bzw. FFT4, FFT5)nachgeschaltet ist, denen Aus- 
gleichsglieder (7, 9; 8, 10; 11; 17, 18) fur den Pegel-. 

20 ausgleich zugeordnet sind. 

20. Gerat nach einem der Anspriiche 13 bis 19, da- 
durch gekennzeichnet, daB zur 
Ermittlung des Spektrums (|F a | 2 ) der Gegenkomponente ein 

25 Differenzbildner (26, 27) zur Bildung der Differenz 
zwischen dem Gesamtspektrum (von FFT1 ) und dem bereits 
isolierten Spektrum (jF^I 2 ) vorhanden ist. 

21. Gerat nach den Anspruchen 18 und 20, d a - 

30 durch gekennzeichnet, daB bei 
schrittweiser Analyse aus zwei entgegengesetzten Fre- 
quenzrichtungen wenigstens zwei Differenzbildner 
(26 , 27) vorhanden sind. 

35 22. Gerat nach einem der Anspriiche 13 bis 21, d a - 
durch gekennzeichnet, daB zur 
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Durchfuhrung des Verfahrens nach den Anspruchen 9 und 
10 eine Addierschal-tung (29, 30), vorzugsweise digitale 
Addierschaltung, vorhanden 1st. 

5 23. Gerat nach einem der- Anspriiche 13 bis 22, d a - 
durch gekennzeichnet, daBzur 
Durchfuhrung des Verfahrens nach Anspruch 11 Rausch- 
synthetisierer (in den Blocken 32, 33 der Figur 5. 
bzw. in der Ausbildung gemaB Fig. 6) vorhanden sind 

10 

24. GerSt nach einem der Anspriiche 13 bis 23, d a - 
durch gekennzeichnet, daBzur 
Durchfuhrung des Verfahrens nach Anspruch 12 eine 
Rechenschaltung (31) zur Bildung der mittleren Fre- 
15 quenz (f ) sowie gegebenenfalls der StrSmung vor- 
handen ist. 
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